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I  P: 


ABSTRACT 


Flows  over  non-ablating  hypersonic  vehicles  are  con¬ 
sidered  when  all  electrons  are  produced  by  aerodynamic 
heating  in  the  boundary  layer.  Detailed  calculations  of 
the  ion  and  electron  distributions  are  presented  for  a 
vehicle  with  a  conical  nose  and  a  semivertex  angle  of 
10  deg,  for  Mach  numbers  8  and  10  at  sea  level.  Except  for 
a  small  region  at  the  nose,  the  boundary  layer  is  turbulent. 
Diffusion  of  electrons  with  respect  to  the  ions  and  neutrals 
is  considered  and  ambipolar  diffusion  is  found  to  occur 
across  most  of  the  boundary  layer  only  at  large  ion-electron 
densities. 
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1.  INTRODUCTION 


The  accuracy  with  wh'^h  a  hypersonic  vehicle  .receiving 
microwave  radiation  through  a  dielectric  window  can  track 
a  target  is  affected  by  electrons  produced  by  aerodynamic 
heating  in  the  boundary  layer  of  the  tracking  vehicle. 

The  plasma  distorts  an  incoming  signal  by  changing  its 
propagation  velocity  through  the  plasma,  refracting  it  at 
the  plasna  radonie  interface,  and  causing  additional  changes 
in  the  passage  of  the  signal  within  the  radome  itself. 
Accordinglv,  wherever  the  electron  number  density  exceeds 
lCq/cm3,  the  calculation  of  boresight  error  requires  a 
description  of  the  electron  distribution.1 * 

In  the  present  case,  detailed  calculations  for  the  ion 
and  electron  distributions  are  given  for  a  tracking  vehicle 
having  the  form  of  a  10-degree  semivertex  cone  flying  at 
zero  angle  of  incidence,  sea  level,  and  Mach  numbers  8  and 
10.  We  assume  that  the  vehicle  is  subjected  to  intense 
aerodynamic  heating  for  a  short  time  duration,  so  that  no 
ablation  occurs  and  the  vehicle  shape  is  unaltered.  Hence, 
with  a  pointed  or  slightly  rounded  non-ablating  nose,  the 
electron  production  occurs  entirely  within  a  small  fraction 
of  the  hot  boundary  layer  adjacent  to  the  vehicle  surface. 
For  the  above  conditions,  the  boundary  layer  is  turbulent 
everywhere  except  for  a  very  small  region  near  the  nose. 

The  random -translational  energy  of  the  electrons  is  at 
most  .0.4  electron-volt,  and  those  electrons  are  subject  to 
large  accelerations  by  small  electric  fields.  Thus,  the 
electron  production  (as  affected  by  recombination)  as  well 
as  the  distribution  are  sensitive  to  the  small  electric 
fields  produced  by  charge  separation,  caused  by  the  faster 
diffusion  of  electrons  relative  to  ions  from  the  region  of 
higher  ionization  production.  The  electric  field  sets  up 
a  conduction  current  which  tends  toward  ambipolar  diffusion 
with  increasing  levels  of  ionization. 

For  the  above  flows,  the  enthalpy  distribution  is 
decoupled  from  the  ionization,  since  the  ionization  enthalpy 
is  insignificant  in  comparison  with  the  local  enthalpy. 

The  ion  mass  fraction  nowhere  exceeds  10  .  Accordingly, 

all  electrons  are  produced  through  the  reaction3 

1  Private  comrnunicat ions, T.B. A. Senior,  Univ.  of  Michigan, 
and  G.  Tricoles,  General  Dynamics-Electronics, 

San  Diego,  California,  March  1971. 

3  J.Hilsenrath,  and  M. Klein,  Tables  of  Thermodynamic 
Properties  of  Air  in  Chemical  Equilibrium  Including 
Second  Virial  Corrections  from  1500° K  to  15, 000° K, 

Arnold  Engineering  Development  Center,  AEDC-TR-65-58, 

March  1965. 
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N  +  0  -  N0+  +  e 


(1) 


Moreover,  the  effect  of  Coulomb  forces  on  collisions  between 
neutral  particles  and  between  ions  and  neutral  particles  is 
negligible.  Consequently,  the  neutral  gas  viscosity  and 
ion  diffusion  coefficients  are  not  affected  by  the  electric 
field. n 


In  the  present  calculations,  we  account  for  the 
turbulent  boundary  layer  by  assuming  a  power-law  velocity 
distribution  for  the  neutral  gas  and  an  associated  iso- 
energetic  enthalpy  boundary  layer.  This  procedure  allows 
for  the  complete  determination  of  the  neutral  gas  flow. 

Binary  and  turbulent  eddy  diffusion  coefficients  are  in¬ 
cluded  in  the  calculation  of  the  ion-electron  drift  velo¬ 
cities.  The  boundary  layer  equations  governing  the  electron 
and  ion  distributions  are  then  solved  exactly  for  dissocia¬ 
tion  equilibrium  and  ionization  non- equilibrium. 

The  uncertainties  in  the  boundary  conditions,  such 
as  body  temperature  and  body  electric  potential,  justify 
use  of  this  approximate  method  of  solution.  These  con¬ 
ditions  are  usually  known  only  to  within  a  range  of  values 
that  can  yield  variations  in  ion  and  electron  distributions 
easily  exceeding  the  error  introduced  by  our  method  of 
solution.  In  addition,  the  effect  on  the  ion  and  electron 
distributions  by  the  neutral  particle  velocity  distribution 
and  boundary  layer  thickness  can  be  readily  obtained  by  our 
procedure.  For  prescribed  boundary  conditions,  the  present 
procedure  permits  a  readily  obtainable  approximate  evaluation 
of  the  plasma  distribution  for  a  wide  range  of  body 
geometries  and  flow  conditions. 

2.  NEUTRAL  GAS  FLOW  AND  THERMODYNAMIC  PROPERTIES 


For  the  cone  of  semivertex  angle  10  deg  considered 
here,  significant  shock  wave  ionization  does  not  occur. 
High  temperature  and  associated  ionization,  say  lOs/cm8 
or  larger,  occur  only  within  a  small  fraction  of  the 
boundary  layer.  As  previously  stated,  all  ions  and 
electrons  are  produced  by  the  reaction  given  by  equation 
(1) .  Assuming  ion  formation  from  Ns  and  08 ,  approximately 
10  ev  per  ion  are  required.  The  energy  required  to  form 
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I.Pollin,  The  Stagnation-Point  Langmuir  Probe  in  a 
Shock  Tube-Theory  and  Measurements,  Harry  Diamond  Labs, 
HDL  TR-1103,  28  June  1963. 


1014/cma  electrons,  which  exceeds  the  equilibrium  ion 
density  formed  on  the  above  body  at  Mach  number  10  and  sea 
level,  is  1.6  (10—* )  joule/cm3.  This  energy  is  insig¬ 
nificant  in  comparison  with  the  Mach  number  10  sea  level 
free  stream  total  enthalpy  of  7.5  joule/cm3.  Thus,  the 
neutral  gas  enthalpy  is  decoupled  from  the  ionization. 


We  assume  the  boundary  layer  neutral  gas  velocity 
profile  corresponding  to  turbulent  flow  to  be  of  the 


form4  -e 


u  =  u  T'm 

5 


(2) 


The  extent  of  the  laminar  sublayer  adjacent  to  the 
vehicle  surface  is  about  1  percent  of  the  total  boundary 
layer  thickness.  Except  for  this  displacement,  the 
turbulent  layer  is  unaffected  by  its  presence,  since  only 
a  small  fraction  of  the  momentum  lies  in  the  sublayer. » »8 9 10 
The  velocity  power  law  given  by  equation  (2)  is  represent¬ 
ative  of  a  large  class  of  fully  developed  turbulent  boundary 


4  R.K.Lobb,  E.M. Winkler,  and  J.Persh,  Experimental 
Investigation  of  Turbulent  Boundary  Layers  in  Hypersonic 
Fiow,  J. Aeronaut .Sci. 22,  1-9-50  (1955). 

5  D.M.Bushnell,  and  I.F.Eeckwith,  Calculation  of  Hypersonic 
Turbulent  Boundary  Layers  and  Comparison  With  Experimental 
Data,  AIAA  Paper  No.  69-684,  June  16-18,  1969. 

6  H.L.Dryden,  Transition  from  Laminar  to  Turbulent  Flow, 

High  Speed  Aerodynamics  and  Jet  Propulsion,  Vol.Vs 
Turbulent  Flows  and  Heat  Transfer,  Princeton  Univ.  Press, 
1959,  Edited  by  C.C.Lin. 

7  C.M. Jackson,  and  R.S. Smith,  A  Method  for  Determining  the 
Total  Drag  of  h  Pointed  Body  of  Revolution  in  Supersonic 
Flow  with  Turbulent  Boundary  Layer,  National  Aero¬ 
nautics  and  Space,  NASA  TN  D-5046,  March  1969. 

8  D.L.Brott,  W.J.Yanta,  R.L.Voisinet ,  and  R.L.Lee.  An 
Experimental  Investigation  of  the  Compressible  Turbulent 
Boundary  Layer  with  a  Favorable  Pressure  Gradient, 

AIAA  Paper  No.  69-685,  June  16-18,1969. 

9  Ibid. 

W.D. Hayes  and  R.F .Probstein,  Hypersonic  Flow  Theory, 
pp  327-339,  Academic  Press, Inc.,  1°59. 
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layer  flows  at  high  and  low  Mach  number,  and  therefore 
provides  an  adequate  basis  for  determining  the  neutral  gas 
velocity.  The  appropriate  value  of  m  appears  to  be  between 
6  and  8;  m  -  8  was  used  in  the  present  calculations9- 1 2 •  As 
previously  stated,  the  effect  on  the  electron  distribution 
caused  by  variation  of  m  within  this  range  is  small  in 
comparison  with  that  produced  by  other  uncertainties,  such 
as  the  value  of  the  vehicle  temperature. 

The  boundary  layer  thickness  5=6 (x) (defined  as  the 
distance  measured  normal  to  the  wall  where  the  velocity  is 
0.99  of  the  local  free  stream  value)  depends  on  such 
factors  as  the  wall  temperature,  vehicle  contour  and  flight 
conditions.  i°-i*  For  example,  at  high  Mach  numbers,  the 
effect  of  insulating  the  body  is  an  increase  in  the  boundary 
layer  temperature  and  thereby,  a  reduction  in  the  density 
and  an  increase  of  the  boundary  layer  thickness.  Usually, 6 
increases  as  the  flow  proceeds  downstream,  and  thereby  new 
fluid  continually  enters  the  boundary  layer  and  is  heated. 

For  fully  developed  turbulent  boundary  layer  flow, 
at  Mach  numbers  1.8  and  4.25  the  boundary  layer  thickness 
on  an  axisymmetric  cone  with  the  semivertex  angle  of 
12-1/2  deg  was  everywhere  found  to  be  less  than  10  percent 
of  the  local  cone  radius  at  all  angles  of  incidence  up  to 
at  least  15  deg.  At  the  latter  angle,  the  value  of  6  on 
the  leeward  side  was  found  to  be  a  factor  of  5  larger  than 
that  on  the  windward  side  at  the  same  downstream  station.13 
The  large  variation  of  boundary  layer  thickness  is  brought 
about  primarily  by  the  circumferential  pressure  gradient. 

A  favorable  pressure  gradient  (pressure  decreasing  in  the 
direction  of  the  flow)  always  tends  to  reduce  the  rate  of 
increase  in  5 ;  the  opposite  is  also  true.  Thus,  at  a 
given  axial  station,  the  boundary  layer  thickness  can 
circumferentially  vary  by  a  factor  of  2  at  2  deg  angle  of 
incidence  and  by  a  factor  more  than  3  at  5  deg.13  Con¬ 
sequently,  unless  the  flow  conditions  and  body  geometry 
are  known,  only  approximate  values  of  6  =6 (x)  can  be 
prescribed. 

11  R.K.Lobb,  E.M. Winkler,  and  J.Persh,  Experimental 

Investigation  of  Turbulent  Boundary  Layers  in  Hypersonic 
Flow,  J.Aeronau.Sci.  22,  1-9-50,1955. 

13  D.M.Bushnell,  and  I .F. Beckwith,  Calculation  of  Hypersonic 
Turbulent  Boundary  Layers  and  Comparison  with  Experimental 
Data,  AIAA  Paper  No.  69-684,  June  16-18,  1969. 

13  W. J.Rainbird,  Turbulent  Boundary  Layer  Growth  and  Sepa¬ 
ration  on  a  Yawed  12  1/2°  Cone  at  Mach  Numbers  1.8 
and  4.25.  AIAA  Paper  No.  68-98,  Jan  22-24,  1968. 
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The  peak  boundary  layer  electron  density  is  reduced 
by  diffusion,  moreover,  the  ionization  process  is  com¬ 
paratively  slow,  so  that  equilibrium  values  of  peak  electron 
density  are  nowhere  attained.  Therefore,  in  an  iso- 
energetic  boundary  layer,  this  means  that  the  peak 
electron  density  will  tend  to  increase  with  increasing  5 . 

For  moderate  angles  of  incidence  a  reasonable  estimate  of 
6  is  given  by  flat  plate  theory  and  Rainbird's  data,1- 
from  which  we  deduce 

6  =  bxa  (3) 

with  a  =  0.8.  The  effect  of  f  on  the  electron  distribution 
can  be  calculated  by  adjusting  b  over  the  range 
0,02  •>  b  >  0.002.  The  smaller  value  of  6  is  associated 
with  low  altitudes,  the  windward  side  of  the  vehicle,  etc. 

-In  the  calculations,  b  =  0.02  was  used. 

Th  pressure  distribution  across  the  boundary  layer 
is  given  by  3p/^y  =0.  For  cones,  in  the  inviscid  region 
p  =  constant  along  a  conical  generator.  Hence,  p  = 
constant  throughout  the  boundary  layer  for  conical 
vehicles.  Inside  the  boundary  layer, 

p  =  p/RTZ  (4) 


For  a  fully  developed  turbulent  flow  in  dissociative 
equilibrium  on  a  conical  body,  the  boundary  layer  enthalpy 
distribution  is  given  by  the  Crocco  relation 


u 


H  =  h  +  (H  -  h,  )  - 
u  fi  13'  u 


when  the  laminar  and  turbulent  Prandtl  and  Lewis  numbers  are 
unity.  i4,ls  For  H  =  h.  ,  we  obtain  an  isoenergetic 
boundary  layer,  so  'that 


(5) 


Throughout  the  boundary  layer,  v  «  u ;  consequently,  the 
v  term  is  not  included  in  the  calculation  for  h. 


14  W.D. Hayes  and  R.F.Probstein,  Hypersonic  Flow  Theory, 
pp  327-339,  Academic  Press  Inc.,  1959. 

55  L.H.Back  and  R.F.Cuffel,  Relationship  Between  Temperature 
and  Velocity  Profiles  In  A  Turbulent  Boundary  Layer 
Along  A  Supersonic  Nozzle  With  Heat  Transfer. 

AIAA  Journal,  Nov„,  1970. 
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Actually,  the  equilibrium  value  of  hb  is  limited  by 
the  recovery  factor,  which  has  the  value  0.90  for  turbulent 
boundary  layers.  The  value  of  h  given  by  (5)  is  every¬ 
where  at  most  10  percent  larger  than  that  obtained  by 
assuming  =  .90  h(  m  the  Crocco  relation.  The  resulting 
difference  in  election  distribution  is  well  within  the 
difference  corresponding  to  the  uncertainty  in  describing 
actual  flight  conditions.  The  effect  of  assuming  an 
isoenergetic  boundary  layer  is  ec[uivalent  to  attaining 
a  given  peak  electron  density  at  a  slightly  lower  Mach 
number  with  a "slightly  different  electron  distribution. 


With  the  .use  of  equation 
becomes 


In  equation  (6),  the  recovery 
.  m 

-  =  (0.1 )F  The  enthalpy  then 

. 01 

(0.1)**  £  n>0. 


(2),  the  energy  equation 

(6) 

enthalpy  is  reached  at 
was  assumed  constant  for 


The  flow  thermodynamic  properties  at  y  =  6  were 
calculated  for  the  conical  flow  at  zero  angle  of  incidence 
with  the  help  of  the  Sims*  tables.1 3 The  formation  of 
equilibrium  values  of  dissociated  N  and  0  takes  place 
within  several  hundred  collisions.17,18,19  The  largest 
mean  free  path  in  the  present  calculations  is  of  the 
order  lO^cm.  Thus,  the  boundary  layer  will  be  in  chemical 
equilibrium  with  respect  to  the  dissociated  N  and  0. 

Since  T  =  T (■" )  and  p  =  constant,  the  N  and  0  distributions 
are  functions  of  ^  only. 


10  J.L.Sims,  Tables  for  Supersonic  Flow  Around  Right 
Circular  Cones  at  Zero  Angle  of  Attack,  National 
Aeronautics  and  Space  Administration,  NASA  SP-3004, 1964. 

1  M.H. Bloom  and  M.H. Steiger,  Inviscid  Flow  with  Non¬ 
equilibrium  Molecular  Dissociation  for  Pressure  Distri¬ 
butions  Encountered  in  Hypersonic  Flight,  Journal  of 
the  Aero/Space  Sciences,  Vol.27,  No. 11,  November  1960. 

10  K.L.Wray,  Chemical  Kinetics  of  High  Temperature  Air 
International  Hypersonics  Conference,  M.I.T., 

Cambridge,  Mass.,  Aug.  16-18,1961. 

19  S.C.Lin,  and  J.D.Teate,  Rate  of  Ionization  Behind  Shock 
Waves  in  Air,  II.  Theoretical  Interpretation,  Avco 
Research  Report  115,  Sept. 1962. 
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Moreover,  because  the  ionization  enthalpy  is  small 
in  comparison  with  the  free  stream  enthalpy,  we  may  assume 
dissociative  equilibrium  in  the  boundary  layer,  and  thereby 
h  *  h(p,T) .  Then,  with  the  use  of  equations  '  (4)  and  (6)  , 
the  thermodynamic  properties  were  found  for  y  <  5  using 
the  tables  of  Hilsenrath  and  Klein,.80  In  every  case,  T, 
is  so  small  that  dissociation  is  zero  and  Z  =  1  at 
y  =  6. Within  the  boundary  layer,  Z  **  Z(p,  T) . 


The  component  of  the  neutral  gas  velocity  normal  to  the 
conical  surface,  v,  is  obtained  from  the  mass  conservation 
equation  for  axi symmetrical  flow, 

(pur)x+(Dvr)y  =  0  (7' 

Since  p  and  u  are  known  f  met  ions  of  T|,and  r  and  £  are 
known  functions  of  x,  the  solution  of  equation  (7)  becomes 


m-H 

au  5T|  “m" 


-  [l+a] 


6U  P 

A 


"T)  *>.  m 


Since  r«x  and  6«x,  by  equations  (2)  and  (8),  v«u; 
therefore,  as  previously  noted,  the  effect  of  v  on  h  is 
everywhere  negligible.  However,  v  has  significance  in  the 
computation  of  the  ion  and  electron  velocities  and  thereby 
on  their  corresponding  density  distributions. 

3.  GOVERNING  EQUATIONS  FOR  THE  CHARGED  PARTICLES 

For  axisymmetric  flow,  the  conservation  equations  for 
the  ions  and  electrons  are 


pusv  +  p  vs. 


(pDiy  +  PDTSy)y  +  (°Ks<sy)y  + 


puix  +  pvsy  =  (pDsy  +  pDTsy)y  -  (pKs<jy)y+  -r 

where  c  ~  c  =  c  and 
e  e  e 


X  -  -  -  lx  =  K  <N>  <0>  -  K  in, 

a  m  dt  1  r 


20  J. Hilsenrath  and  M. Klein,  Tables  of  Thermodynamic 
Properties  of  Air  in  Chemical  Equilibrium  Including 
Second  Virial  Corrections  from  1500°  K  to  15.,  000°  K, 
Arnold  Engineering  Development  Center,  AEDC-TR-65-58, 
March  1965. 
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and  the  formation  and  recombination  rate  constants  are 
given  as"1 #a2 »23 >s* 


Kf  I:  5(10  exp (-32, 50^/T)  and 

(12) 

K  =  3(10-3)T-1*5  cm2 /ion-sec, 
r 


Fig.  1(a)  Boundary  layer  parameters 
D,  P  and  T  -  mach  number  8 


Fig.+  l(b)  Boundary  layer  parameters 
D,  P  and  T  -  mach  number  lu 


Sl  S.C.Lin,  and  J.D.Teare,  Rate  of  Ionization  Behind  Shock 
Waves  in  Air,  II,  Theoretical  Interpretation,  Avco 
Research  Report  115,  Sept. 1962. 

22  M.A.Biondi,  Recombination  Processes  (Charged  Particles), 
Defense  Atomic  Support  Agency  Reaction  Rate  Handbook, 

DASA  1948,  Oct. 1967,  pp.  11-8  to  11-10. 

23  A.Q.Eschenroeder,  and  T.Chen,  Near-Wake  Ionization  Behind 
A  Sphere  in  Hypersonic  Flight.  I. React ion  Kinetics, 
General  Motors  Defense  Research  Labs.TR  65-01  H. 

Sept. ,1005. 

A.Q.Eschenroeder,  and  T.Chen,  Near-Wake  Ionization  Behind 
A  Sphere  in  Hypersonic  Flight.  II.  Influence  of  Flight 
Conditions,  General  Motors  Defense  Research  Labs., 

TR  65-01  Q.  Dec. 1965. 
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The  diffusion  coefficient  of  NO*  in  air, denoted  by 
D,  is  about  the  same  as  that  for  NO  in  air35  .  Consequently, 
using  the  Leonard-done s  potential,  for  constant  pressure  D 
is  a  function  of  temperature  only.  Since  T  =  T('),  using 
the  values  given  in  Table  V  of  reference  26 ,  the  dependence 
of  D  on  n  is  shown  in  figure  1  for  sea  level  flight  at  Mach 
numbers  8  and  10,  where  the  boundary  layer  pressures  for  the 
two  cases  are  4.29  and  6.00  atmospheres,  respectively. 

The  binary  diffusion  coefficient  for  the  electrons 
in  air  can  be  found  from  the  approximate  relation 

6  =  (iS/M)^  0?/$)  D  =  2j4D,  (13) 


where  we  have  assumed  T  =  T, 


When  the  turbulent  Prandcl  and  Lewis  numbers  are  unity, 
the  turbulent  eddy  kinematic  viscosity  and  turbulent  eddy 
diffusion  coefficients  become  equal2 7  The  mixing  length 
theory  of  Prandtl28,  gives  .$  =  Xy  and  u  'v  (™) "  ,  so 

5  .  With  X  =  0.4  and  (2),  there  results  for 

dy 

both  ions  and  electrons. 


that  Dt  =  t 


Dt  =  0.02  u^  5 


t8/" 


(14) 


Since  uf  3  (10s  )  cm/sec,  for  x>  1cm, 

DT>l20T)r  cm2 /sec.  By  fig.  1,  except  in  the 
neighborhood  of  f]  =  0  (that  is,  except  in  the  laminar 
sublayer),  DT»D.S imilarly,  DT»D  at  somewhat  larger  X. 

For  small  values  of  the  ratio  of  the  electric  field 
to  the  total  gas  pressure,  the  Einstein  equation  K  =  eD/kT  or 

K  =  11,600  D/T  (15) 

25  I.Amdur  and  E. A. Nason,  Properties  of  Gases  at  Very  High 
Temperatures,  Physics  of  Fluids,  Vol.l,  No. 5, 
pp.  370-83,  Sept-Oct  1958. 

36  I.Pollin,  The  Stagnation-Point  Langmuir  Probe  in  a 
Shock  Tube — Theory  and  Measurements,  Harry  Diamond  Labs, 
HDL  TR-1103 ,  28  June  1963. 

37  W.D. Hayes,  and  R.F.Rrocstein,  Hypersonic  Flew  Theory, 
pp  327-339,  Academic  Press  Inc.,  1959. 

38  H.Schlichting-,  Boundary  Layer  Theory, 

McGraw-Hill  Book  Co,  4th  Ed.,  1960,  pp  477-490. 
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relates  the  mobility  and  binary  diffusion  coefficients  for 
both  electrons  and  ions,  ref.89 


The  y  component  of  the  ion  or  electron  current 
density  is  defined  by  the  equation 

J  =+vDen, 


(16) 


(j>0  is  in  the  direction  toward  the  body  surface)  where  the 
y  component  of  the  average  drift  velocity  is 

+  ± 


vd=  _v  -f  DS-/i3  +  Dms../s  ±  Ktf 


y 


T  y 


(17) 


and  where  (+)  is  for  ions  and  (-)  is  for  electrons.  Sub¬ 
stitution  of  equations  (16)  and  (17)  into  equations (9)  and 
/10)  and  w  =  y/f  gives 


fen  us  Wejn 

«  e  x  e  ,  v, 

j  =  -  •  ■  ■  -  —  •  -  en  s  ( — — 

?  c  c  e  P  T 

G  G 


(18) 


for  either  ions  or  electrons. 

Eliminating  vD  in  equations  (16)  and  (17)  gives 


+ 

s 


P  +  ± 
T  B 


+ 

JPf 

-en 


—  T 

SVf 


* 


(19) 


where  the  (-)  and  (+)  signs  preceding  the  electric  mobility 
term  are  for  ions  and  electrons,  respectively. 


Assuming 


?v 


electric  field  is 


->  *>  c  2  * 


10" 


<T1 


8.85 


xx 


en  5' 
e 


,  Poisson's  equation  for  the 


(s  -  s) 


(20) 


Previously,  u,  T  and  P  were  made  functions  of 
As  we  have  seen,  equilibrium  dissociation  may  be  assumed 
and  thereby  the  <N>  and  <0>  are  only  functions  of  t|. 

Since  the  pressure  is  constant  throughout  the  boundary  layer, 
p  and  £  are  also  functions  of  T]  only.  By  (14),  at  a  given 
x,  D  =  D  (ti )  .  The  five  equations  (18) -(20)  are  explicit 
expressions  fcr  the  derivatives  jr  *  ,  it  )_  in  terms  of 

V 1  11  ^ 

J  '  s  '  s  '  C  k  Consequently,  the  computation 


A.  Von  Engle,  Ionized  Gases,  Oxford  Univ.  Press, 1955. 
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for  J,  3  and  $r  varies  with  r  for  different  x  through  v=v(x,^), 

t=c:  (x)  and  %  =|  (x,t<). 

X  X 

4.  B0UND7iRY  CONDITIONS 

The  solution  of  equations  (18), (19)  and  (20)  requires 
specification  of  five  boundary  conditions  in  addition  to 
initial  conditions  at  x  =  0.  For  sharp  nose  cones 

at  x  -  0,  s  =  0  (21) 

As  previously  stated,  all  ionization  production  occurs 
within  a  small  fraction  of  the  boundary  layer  adjacent  to 
the  vehicle  surface.  Whether  a  wall  is  catalytic  or  non- 
catalytic  with  respect  to  charged  particle  recombination 
depends  on  the  gas-solid  surface  interaction,  which  in  turn 
depends  somewhat  on  the  wall  temperature.  For  the  wall 
temperature  of  1000"  K  assumed  by  Wang3-,  experimental 
observation^  agree  with  the  calculations  obtained  using  the 
conditions  s(0)  =  0.  Similarly,  we  assume 

at  *  =  0,  S  =  0  (22) 

The  temperatures  at  the  edge  of  the  boundary  layer  and 
for  some  distance  therein  are  very  low  (see  fig.  1) .  Except 
for  the  diffusion,  this  condition  would  give  s  =  0  at  T>  =  1. 
Calculations  show  that  nonzero  boundary  conditions  for 
S  at  -  =  1  wherein  s_-0  at  v  =  1  principally  affect  only  the 
*  distributions  near  ?  =  i .  Consequently,  we  assume 

at  -  =  1,  s  =  0  (23) 

Tfce  boundary  condition  for  an  insulated  surface  is 
J  t=  j  +  J  =  0,  so  that  (1)  is  a  function  of  x.  Including 
the  passage  of  current  into  the  wake,  ip  (1)  is  not  known  for 
a  conducting  surface.  The  value  of^jb  near  Tj  =  1  is 
necessarily  small  because  of  the  small  values  of  s  near 
^  =  1  and  the  electric  field  would  appear  to  be  small  at 
n  =  1.  Accordingly,  for  all  x,  we  assume 

at  -n  =  1,  0^  «  0  (24) 

Thus,  (24)  is  an  approximate  boundary  condition  for  both 
insulated  and  conducting  surfaces. 


a'o  k.  Wang,  Electron  and  Ion  Distributions  in  Chemical 
Nonequilibrium  Boundary-Layer  Flows,  AIAA  Journal, 
p.316,  April  1969. 
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In  the  next  section,  the  Mach  number  8  calculations 
show+that  is  not  significant  in  the  determination  of  s 
and  J,  since  the  terms  are  small  compared  to  the  other 
diffusion  terms.  ^Tor  Mach  10,  4>n  is  small  for  small  x, 

•and  again  the  effect  of^,  on  the  calculations  is  negligible. 
At  larger  x,  the  faster  electron  diffusion  results  in  signifi¬ 
cant  values  of  at  n  =  0,  fig.  3  {c) .  Upon  using  these 
values  of  for  the  boundary  condition  (24) ,  the  effect  on 
the  s  and  j  distributions  (not  shown)  was  found  to  be 
limited  to  the  neighborhood  of  ij  =  1  and  did  not  affect  the 
rat  -  =  0  and  the  peak 

Finally,  ‘in  the  evaluation  of  0,  we  assume 

at  =  1,  $  =  0  (25) 

5.  NUMERICAL  CALCULATIONS  FOR  THE  ION  AND  ELECTRON  FLOW  IN 
THE  BOUNDARY  LAYER 

The  \  (x,-)  in  (18)  can  be  approximated  by  an  appropriate 
finite  difference  expression  involving  I(x,^),  I (x-^x,m), 
etc.,  so  that  the  system  (18) -(20)  with  the  boundary  con¬ 
ditions  (21) -(25)  ultimately  reduces  to  solving  a  two-point, 
nonlinear  boundary- value  problem  at  each  x.  This  was  solved 
numerically  both  by  integrating  the  differential  equations 
and  by  matrix  methods.31As  determined  by  reducing  the  step- 
size  cf  the  calculation  intervals  /‘x  and  At  at  any  station 
x,  the  accuracy  of  any  quantity  was  found  to  be  within 
about  5  percent  providing  the  quantity  was  not  reduced  to 
less  than  1  percent  of  its  maximum  value.  Improvement  of 
the  calculation  precision  requires  longer  computer  running 
time.31 


The  numerical  data  is  summarized  in  tables  I  and  II  and 
figures  2  and  3.  In  particular,  the  tables  give  the  printout 
of  the  independent  variables  at  typical  stations  0<x^j-oo 
as  well  as  the  contributing  terms  to  v^.  The  graphs  of  r 
were  obtained  from  the  relation 

n  =  nes/p  (26) 

where  n  =6.54  (10l0)  and  5.63  (101?)  for  Mach  numbers 
8  and  e10  and  the  corresponding  P  =  P^T))  functions  are 
given  in  figure  1. 


31  A.  Hausner,  Ion  and  Electron  Distributions  in  the  Boundary 
Layer  of  Hypersonic  Vehicles  For  Chemical  Non-Equilibrium 
Flow  -  Part  II:  Method  of  Solution  and  Computer 
Program,’  Harry  Diamond  Labs.,  HDL-TR  1567  1971. 
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+  + 

All  functions  s,  j,  and  are  zero  at  x  =  0.  The 
enthalpy  and  density  profiles  of  the  neutral  gas  flow  result 
in  the  total  ionization  production  occurring  within  a  small 
fraction  of  the  boundary  layer  in  the  neighborhood  of  the 
vehicle  wall.  For  all  x,  the  diffusion  caused  by  con¬ 
centration  gradient  is  directed  from  -  -  0  toward  n  -  1. 
Thus,  in  the  neighborhood  r  =  0,  ion-electron  recombination 
is  negligible  and  the  ionization  production  remains  con¬ 


stant  for  all  x.^  The  diffusion  effects  become  significant 


with  increasing  s  and  therefore  with  increasing  x. 

because  of  the  diffusion..-^  „  ,  -  .. 

.  ..  .  ,  ^Oandn-'-'n  m  the  ■" 

neighborhood.  „  -  e 


Hence, 
=  0 


s 

l0-« 


/  =0.5  a 

.  ; 

“o  50  100 

X 

Fig.  2(a)  Ion  and  electron  densities 
at  ri  -  0.01,  0.1  and  0.5  mach  number  b 


V 


Fig.  3(a)  Normalized  electric  field  & 
potential  at  X*  100cm  -  mach  number  8 


s 


x 


Fig.  2(b)  Ion  and  electron  densities 
at  "  =  0.01,  0.1  and  0.5  mach  number  10 


0  0  5  1.0 


1 

Fig.  3(b)  Ion  and  electron  densities 
at  X  =  5.62  and  100cm  -  mach  number  8 


Fig.  3(c)  Normalized  electric  field  & 
potential  at  X=5.62  and  100cm  mach 
number  10 


Fig.  3(d)  Ion  and  electron  densities 
at  X  =  5.62  and  100cm  mach  -  roac.h  number  10 


From  previous  work3",  the  u  velocity  is  sufficiently 

small  at  ~  10-4  to  allow  time  for  .ioni  zation  eauilibrium 

'well  before  x  =  100,  particularly  forMach  number  i0. 
Ionization  equilibrium  at  ~  •-  10-4  is  denoted  by  g  =  i. 

By  f igures3 (b) and  3(d),  for  Mach  numbers  8  and  10,  the 
diffusion  reduces  the  peak  electron  density  at  x  =  100  by 
a  factor  of  more  than  10  and  103/  respectively.  On  the 
other  hand,  for  x  =  100,  the  ±  at  "  •>  .1  are  much  larger 

than  those  corresponding  to  ^ionization  equilibrium. 

This  occurrence  is  explained  by  the  diffusion  from  the 
region  '  0.1  and  the  low  recombination  rate. 


For  the  Mach  number  10  calculations  we  observe  the 
following  behavior:  Because  of  the  faster  electron 
diffusion  over  that  of  the  ions,  for  small  x>s  '*) '"<S  (" )  , 
fig.  2(b).  As  x  increases,  s  as  g  near  T  =  1  and  the 
’"-region  of  approximate  ion-electron  equality  progresses 
from  -  1  toward  ”=0.  Finally,  at  x  =  100,  s  s  every¬ 
where  except  for  ~  0.01  (figs. 2 (b)  and  3(d)).  This 

difference  between  i  and  £  gives  rise  to  a  0  which 

increases  sharply  near  n  =  0  and  is  a  monotonicaliy 
increasing  function  of  x  at  ~  =  0  (fig.  3(c)).  The  values 
for  are  obtained  by  integration  of  -’-(fig.  3(c)). 


I.  Poll in,  Plasma  Formation  in  the  Boundary  Layer  of  a 
Sharply  Pointed  Hypersonic  Vehicle,  Harry  Diamond  Labs., 
HDL-TR-i453, June  1969. 
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Referring  to  the  £  development  with  x  for  various 
Tl  in  fig.  2(b),  for  sufficiently  large  x,  steady  state 
values  of  §  will  finally  evolve  for  the  Mach  number  10  case  with 
the  following  characterist ics  : 

Ambipolar  diffusion  occurs  everywhere  except  in  the 
immediate  neighborhood^ f  r  =  0  where  g  <  s*  Because  of 
the  large  the  (s  -  s)  difference  will  give  rise  to  a 
large  m  >  0  at  m  =  0.  The  peak  ion-electron  densities 
occurring  in  the  region  10“4  will  be  about  two  orders 
of  magnitude  below  ionization  equilibrium,  whereas  the 
ionization  for  i>>0.1  will  greatly  exceed  equilibrium 
conditions. 

The  diffusion  coefficients  for  Mach  number  8  are 
slightly  smaller  than  for  Mach  number  10.  For  Mach  number 
8,  the  diffusion  is  more  effective  in  reducing  the  peak  n 
because  of  the  slower  rate  of  ionization  production.  Now, 
the  n  are  everywhere  markedly  reduced  so  that  the  contri¬ 
bution  to  vD  by  the  electric  field  resulting  from  charge 
separation  is  everywhere  negligible e (fig. 3  (a;  and  table  I). 
Hence,  for  Mach  number  8,  the  electron  diffusion  everywhere 
exceeds  that  corresponding  to  ambipoiar  diffusion  and,  for 
all  x,  n  >  lo  n  (fig.2(a)).  (if  the  electric  field  arising 
from  the  charge  separation  was  not  negligible,  a  counter 
flow  diffusion  would  occur  which  would  increase  the  peak  n 
and  reduce  the  peak  *  These  effects,  together  with  the 
lower  n  (by  two  orders  of  magnitude)  result  n  reduced 
peak  n  0f  about  5  to  6  orde  :s  of  magnitude  lower  than  those 
for  Mach  number  10. 


6.  CONCLUSIONS 


The  ion-electron  distribution  in  the  turbulent  boundary 
layer  has  been  calculated  for  a  hypersonic  vehicle  consisting 
of  a  10  degree  semivertex,  sharply  pointed  cone  at  Mach 
numbers  8  and  10  for  sea  level  flight.  Binary  and  turbulent 
eddy  diffusions  across  the  boundary  layer  normal  to  the 
vehicle  surface  act  to  reduce  the  peak  whereas  the  electric 
field  arising  from  the  charge  separation  acts  to  cause  ambi*- 
polar  diffusion  an<|  thereby  acts  to  increase  the  peak  n  and 
decrease  the  peak  n.  For  the  cases  considered,  a  balance 
between  the  ionization  production  and  diffusion  results  in 
the  attainment  of  a  steady  state  peak  electron  density  of 
about  two  orders  of  magnitude  less  than  equilibrium  ionization 
at  Mach  number  10  and  a  reduction  of  about  six  orders  at 
Mach  number  8.  Moreover,  the  distance  along  the  missile 
surface  required  to  attain  the  steady  state  peak  electron 
density  in  the  present  calculation  is  several  times  larger 
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than  the  corresponding  distance  to  attain  peak  equilibrium 
ionization  for  the  same  flight  wherein  all  diffusion  is 
neglected . 

The  calculation  procedure  allows  for  an  easy  evaluation 
of  the  effect  of  the  neutral  gas  enthalpy  profile,  boundary 
layer  thickness,  binary  and  eddy  turbulent  diffusion  co¬ 
efficients,  etc.,  on  the  ion-electron  turbulent  boundary 
layer  distributions. 

The  described  procedure  can  be  used  to  determine  the 
ionization  distributions  in  the  turbulent  boundary  layer 
of  sharply  pointed  vehicles  for  Mach  numbers  up  to  15  at 
sea  level  and  20  at  the  altitude  10s  ft.  The  ion  mass  frac¬ 
tion  exceeds  lO-1*  at  higher  Mach  numbers,  and  the  electron 
production  is  no  lonoer  adeauatelv  aiven  hv  reaction  (1). 

The  assumption  of  the  boundary  conditions  s  =  0  at 
x  =  0  requires  a  sharply  pointed  nose.  For  slightly  rounded 
and  blunt  nosed  vehicles,  in  the  nose  shock  region,  shock 
induced  ionization  governs  the  ion  and  electron  distribu¬ 
tions;  however,  diffusion,  initially  aided  by  recombination, 
will  produce  rapid  reductions  of  the  ^  as  the  flow  pro¬ 
ceeds  downstream.  Consequently,  for  the  extension  to 
blunt  nose  bodies,  we  can  assume  the  initial  conditions 
to  be  those  corresponding  to  equilibrium  ionization  behind 
the  nose  shock  and  proceed  therefrom  as  before. 


7.  LIST  OF  SYMBOLS 


Cj  mass  fraction  of  component  j 

D  binary  diffusion  coefficient,  cms/sec 

Dt  turbulent  eddy  difiusion  coefficient,  cm2 /sec 

e  electron  charge  =  1.60  (10  "  )  coulomb  =  4.80  (lO-50)  esu 

h  enthalpy  per  unit  mass  of  mixture 

H  total  enthalpy  per  unit  mass  of  mixture  =  hv-V3 
J  y  component  of  ion  current  density,  Amp/cm2 
K  electric  field  mobility  coefficient,  cm2/V-sec 
Kf  rate  constant  for  ion-electron  formation,  cm-5 /ion-sec 

Kr  rate  constant  for  ion-electron  recombination,  cm3 /ion-sec 

m  velocity  profile  power  law  coefficient  =  8 
M  mass,  g 

n  ion  or  electron  number  density,  cm  3 
<N>  nitrogen  number  density,  atom/cm3 
<0>  oxygen  number  density,  atom/cm3 
p  pressure  xg/m3  or  atmosphere 
P  density  ratio  =  o&/p 

r  radial  distance  of  vehicle  surface  from  axis  of  symmetry, cm 
R  gas  constant  for  air  =  1.987  cal/molec K 
s  mass  fraction  for  ions  or  electrons  =  c/c 

e 

t  time,  sec 

T  absolute  temperature,  °K 

u  x  component  of  neutral  gas  mixture  velocity,  cm/sec 

v  y  component  of  neutial  gas  mixture  velocity,  cm/sec 

V  speed  of  neutral  gas  mixture  velocity  =  /u3  +  v2 ,  cm/sec 
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W  mass  rate  of  formation  of  ions  or  electrons,  9/300-011^ 
x  distance  along  meridian  profile,  cm 
y  distance  normal  to  the  surface,  cm 
Z  compressibility  factor  =  p/oRT  (Z  =  1  at  STP) 
f  boundary  layer  thickness,  cm 

dimensionless  boundary  layer  thickness  -  y/5 
k  Boltzmann  constant  ~  1.38  (10-3“ ) joule/0 K 

:  mass  density,  kg/m? 

r  electric  potential,  V 


Subscripts 


b  vehicle  surface 

D  y  component  of  total  ion  or  electron  velocity,  cm/sec 
e  recovery  temperature  condition  ^  10-4 ) 
x  partial  differentiation  along  meridian  profile 
y  partial  differentiation  along  surface  normal 

'  edge  of  boundary  layer 

■”  partial  differentiation  along  dimensionless  surface  normal 
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